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ABSTRACT 

We explore the galaxy formation physics governing the low mass end of the HI mass 
function in the local Universe. Specifically, we predict the effects on the HI mass 
function of varying i) the strength of photoionisation feedback and the redshift of 
the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback 
prescription, and iv) the efficiency of star formation. We find that the shape of the 
low-mass end of the HI mass function is most affected by the critical halo mass below 
which galaxy formation is suppressed by photoionisation heating of the intergalactic 
medium. We model the redshift dependence of this critical dark matter halo mass by 
requiring a match to the low-mass end of the HI mass function. The best fitting critical 
dark matter halo mass decreases as redshift increases in this model, corresponding to 
a circular velocity of ~ 50kms“^ at z = 0, ^ 30kms“^ at z ~ 1 and ~ 12kms“^ 
at z = 6. We find that an evolving critical halo mass is required to explain both the 
shape and abundance of galaxies in the HI mass function below Mhi lO®ft.“^M0. 
The model makes specific predictions for the clustering strength of Hl-selected galaxies 
with HI masses > and > IO^/i“^M0 and for the relation between the HI 

and stellar mass contents of galaxies which will be testable with upcoming surveys with 
the Square Kilometre Array and its pathfinders. We conclude that measurements of 
the HI mass function at z ^ 0 will lead to an improvement in our understanding of 
the net effect of photoionisation feedback on galaxy formation and evolution. 

Key words: galaxies: formation - evolution - large-scale structure of the Universe 
- radio lines: galaxies 


1 INTRODUCTION 

There is a wealth of observational evidence that the Uni¬ 
verse experienced an “Epoch of Reionization” (EoR) at 
z ^ 6 (e.g. Ouchi et al. 2010 Mesinger et al. 2011| |Mc- 


Greer et al.| 2011 Shull et al.||26l2 | dnring which time the 

cosmic diffuse neutral hydrogen was “re-ionized” by a back¬ 
ground of ultra-violet (UV) and X-ray radiation produced 


by the first generation of stars and galaxies (e.g. Barkana 


& Loeb 2007 Robertson et al. 20101. This photo-ionizing 


background heated the gas in low-mass dark matter halos, 
leading to its expulsion from their shallow potential wells 
( Shapiro et al.pOOd l. In addition, an increased Jeans mass 
prevented infall of gas onto low mass halos with virial tem¬ 


peratures 10'* K (cf. Efstathiou 19921. This implies that 


there should be a critical dark matter halo mass - and there¬ 
fore a halo circular velocity - below which galaxy formation 
should be suppressed during reionization. Recently, |Brown| 
et al. (20141 presented high-resolution spectroscopy of stars 


in six ultra-faint dwarfs within the vicinity of the Milky- 
Way and analysed their stellar population properties. Brown 
et al. concluded that these galaxies formed ~ 80% of their 
stars by z ~ 6, and therefore are the best nearby candidates 
for probing star formation quenching during the EoR. These 
galaxies have maximum velocities ~ lOkms”*, imposing an 
observational limit for star formation quenching due to the 
reionization. 


hansikk@unimelb.edu. au 


Rees (19861 used analytical arguments to deduce that 
this critical circular velocity should be Uut — 30kms“*, 
reasoning that the rate at which gas is heated by an ion¬ 
izing background is balanced by the rate at which it cools 
radiatively. Subsequent studies have tried to sharpen the 
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2 Han-Seek Kim et al. 


prediction for Vcut, with calculations based on semi-analytic 
approaches (e.g. Efstathiou 1992 Benson et al.||2002l, ide¬ 


alised ID hydrodynamical calculations (e.g. Thoul & Wein¬ 


berg 19961 and fully cosmological hydrodynamical simnla- 
tions (e.g. Gnedin 2000| Okamoto et al. 20081. However, 


there is still debate about the precise value of ■ Recently, 


Sobacchi & Mesinger (20131 used ID collapse simulations to 


show that the critical halo mass depends on the redshift. 


following earlier work by Dijkstra et al. (20041. 


The HI mass function in the local Universe provides in¬ 
formation about galaxy formation physics in low mass dark 


matter halos. Indeed, Kim et al. (20131 recently showed that 


feedback effects are important in determining the HI mass 
function in the local Universe. In particular, the low mass 
end of the HI mass function is governed by photoionisa¬ 
tion feedback processes in small dark matter halos. In this 
paper we discuss the question of whether such modelling 
can be used to gain further insight into the value of Ucut- 
Specifically we explore the effect of photoionisation feed¬ 
back on the HI mass function using the GALFORM semi- 


analytic galaxy formation model (Cole et al. 20001. This 


model reproduces the observed HI mass function at z=0, 
accurately matching its amplitude and shape at interme¬ 
diate and high HI masses. We compare the models with 
observational estimates of the HI mass function from the 
HIPASS (HI Parkes All-Sky Survey; see [Meyer etH][2004t 
and ALFALFA (Arecibo Legacy Fast ALFA Survey; see jGio-j 
vanelli et al.|200'5l surveys. We use the version of the semi- 


analytical galaxy formation model GALFORM developed by 
[Lagos et al!] ( [20TTj [Ml] ) as the default model. 

The structure of the paper is as follows. In © we pro- 
vide a brief overview of the Lagos et al. (2011 \ model. In !|^ 
we show how the different ingredients influence the HI mass 
function. In !|^we describe the modified model for photoion¬ 
isation feedback. We finish with a summary of our results in 


2 THE GALAXY FORMATION MODEL 

We use the semi-anal 3 rtical galaxy formation model GAL¬ 
FORM (cf . [Cole et al.|[2000[ |Baugh[[2006[ ) described in [La- 


gos et al. (2011 20121 to predict the properties of galaxies 
forming and evolving in the ACDM cosmologjQ GALFORM 
models the key physical processes of galaxy formation, in¬ 
cluding the gravitationally driven assembly of dark matter 
halos, radiative cooling of gas and its collapse to form cen- 
trifugally supported discs, star formation, and feedback from 
supernovae (SNe) and active galactic nuclei (AGN). 


et al.|2009l with a variety of observations (e.g. the HI mass 


function, ^‘‘CO (1-0) luminosity function, and correlations 
between the ratio H 2 /HI and stellar and cold gas masses in 
Lagos et al. ( 2012[)) a nd found that the empirical law of Blitz 
& Rosolowsky (20061 (see also Leroy et al.|2008' I is favoured 
by these data. This law is of the form 


EsFR = Z^SF fmol Eg 


( 1 ) 


where Esfr and Egas are the surface densities of the star 
formation rate (SFR) and total cold gas mass respectively, 
j^SF is the inverse of the SF timescale for the molecular gas 
and fmol = Emoi/Egas = Rmoi/(Rmoi + 1) is the molecular 
to total gas mass surface density ratio. Rmoi is defined by 
Rmoi = E(iL 2 )/E(H/) = (Pext/.Po)'^'’”'"’- -Pext is the inter¬ 
stellar gas pressure in a galaxy. We calculate Pext using the 
approximation from Elmegreen] ([1993|), 


Pext ~ 2 ^ 


I E* 

a. 


( 2 ) 


where G is the gravitational constant and E* is the to¬ 
tal surface densities of the stars. The velocities dispersions 
of the gas and stars are given, respectively, by Ugas and 
at. Th e va lues allowed by th e observations of [Leroy et al.[ 
(20081 and Bigiel et al. (2010| are i/sf = 0.52 ± 0.25 Gyr“ , 
log(Po/kB[cm-'K]) = 4.54 ±0.07 and = 0.92 ±0.07. 
Importantly for the work we present in this pa per, [Lagos[ 
et al. (20111 showed that the Blitz & Rosolowsky (2006 \ law 


is able to broadly reproduce the HI mass function at z=Q at 
high and intermediate HI masses, Mhi > 10®/i“^M©. This is 
because it suppresses star formation in lower mass galaxies 
due to their low Egas, thereby reducing SNe feedback and 
allowing these galaxies to retain larger gas reservoirs. Note 


that we use the Lagos et al. (20121 model as the default 


model in this paper. In addition to the [Blitz fc Rosolowsky[ 
(20061 star formation law, this model has longer duration 


starbursts compared with the model of Bower et al. (20061. 


As discussed in [Kim et al.[ ( [2013[ ) , photo-ionization is 
predicted to have a dramatic impact on star formation in 
low-mass galaxies. This is because the presence of a photo- 
ionizing background both modifies the net cooling rate of 
gas in halos by removing the “hydrogen peak” in the cool¬ 


ing curve (cf. Fig. 1 of Benson et al. 20021 and increases 
the temperature of the intergalactic medium (IGM) such 
that its thermal pressure prevents gravitational collapse onto 
low-mass halos (e.g. Efstathiou|1992 Okamoto et al.|2008 |. 
As a result, only those halos that hosted cold gas prior to 


Lagos et al. ( 2011 1 extended GALFORM by modelling re-ionization can form stars (e.g. Hoeft et al. 20061. GAL- 


the splitting of cold gas in the interstellar medium (ISM) 
into its HI and H 2 components and by explicitly linking star 
formation to the amount of H 2 present in a galaxy. [Lagos[ 
[et al.[ ( [2011[ ) compared empirically and theoretically derived 
star formation laws (cf. Blitz fc Rosolowsky|2006' Krumholz 


FORM includes the Benson et al. (2002 \ prescription for 


suppressing the cooling of halo gas onto the galaxy - this 
occurs if the host halo’s circular velocity I4irc lies below a 
threshold Kut at redshift Zcut- The values in the default |La-| 


gos et al. (20121 model are Ucut = 30kms and ^cut = 10. 


The default Kut is motivated by the results of hydrody- 


^ Recall that the cosmological parameters adopted for the Mil- 


lennium Simulation are the total matter density Hm = 0.25 
the baryon densi^~*nb~*=~~0±4±~The‘^acuunr*energ^*density 

Ha = 0.75, the Hubble parameter Hq = 73kms“^ Mpc“^, the 


pnmordTar"scaTar'spectranndejrTrr'="^"anmie'^uc±mtToir^nR 

plitude (78 = 0.9. 


namical simulations by Hoeft et al. (20061 and Okamoto 


et al. (20081. However, several studies indicate that I4ut 
should be redshift dependent (Dijkstra et al.|2004 Sobacchi 


fc Mesinger[[2013 1. The implementation of a single value of 


Zcut and Ucut may therefore be overly simplistic. We revise 
this modelling in this paper. 
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The HI mass function as a probe of photoionisation feedback on low mass galaxy formation 3 


3 UNDERSTANDING THE THE LOW MASS 
END OF HI MASS FUNCTION IN THE 
LOCAL UNIVERSE 


The self-consistent treatment of HI in Lagos et al. 1 2011 
| 2012 [ ) results in a model that shows good agreement between 
model predictions and the HIPASS and ALFALFA observa¬ 
tions for HI masses greater than Mhi ~ lO®h“^M 0 . How¬ 
ever the predicted HI mass function between and 

lO®/i“^M 0 (i.e., the low mass end of the HI mass function) 
fails to capture the shape and abundance of the observed HI 
mass function. In this section we investigate how predictions 
of the low mass end of the HI mass function at 2 ~0 are influ¬ 
enced by the assumed photo-ionizing background feedback, 
cosmology, supernovae feedback prescription, and star for¬ 
mation law efficiency. We also introduce a description of the 
possible redshift dependence of the photo-ionizing feedback 
strength. Throughout the paper we use Monte-Carlo gener¬ 
ated merger histories (Parkinson et al. 20081. We adopt a 
minimum halo mass of Mhaio.min = 5 x 10“ h~^ M 0 . From 
the chosen minimum halo mass, we can predict the neutral 
hydrogen gas mass structures down to the current observa¬ 
tional limits (Martin et al. 12010 1 . 


3.1 Photoionisation feedback at the end of 
reionization 

In the standard implementation of GALFORM, the effect 
of photoionisation feedback induced by the epoch of reion¬ 
ization is modelled by imposing a circular velocity cut¬ 
off Vcut on galaxy formation at a redshift of the end of 
reionization 2 cut. In Fig. we show the resulting HI mass 
function assuming this standard model for reionization for 
a range of photoionisation strengths Vut (at hxed 2cut) 
and redshifts marking the end of reionizaiton 2 cut (at fixed 
Kut = 30kms“^). These models can be compared with ob¬ 
servations from ALFALFA ( Martin et al.|2010| and HIPASS 
(Zwaan et al. 20051. The single values of I4ut and 2cut 


introduce a non-monotonic feature in the HI mass func¬ 
tion that is not present in the observations. The top panel 
of Fig. shows that the predicted abundance of galaxies 
with 2 x 10 ^/ 1 “^ M 0 < Mhi < lO®/i“^M 0 is larger than 
observed when we choose the Vut = 25kms“^, which 
matches the abundance at Mhi < 10 ^ h~^ M 0 (red dotted 
line). In contrast, the predicted abundance of galaxies with 
Mhi < 10 ^ h~^ M 0 is lower than observed when we select 
Vcut = 35kms“^, which fits the abundance of galaxies with 
2 X lO^/i“^M 0 < Mhi < 10® h~^ M 0 . This is evidence of the 
simple modelling of photoionisation feedback leading to the 
non-monotonic feature discussed above. 

We also investigate the effect of varying the redshift of 
reionization at a fixed Mut = 30kms“^. The bottom panel 
of Fig. shows that the value of 2 cut does not influence 
the low mass end of the HI mass function for values rang¬ 
ing between 2cut = l to 10 (in the case of Kut = 30kms“^) 
indicating that the HI mass function in the local Universe 
cannot be used as a probe of the end of the epoch of reioniza¬ 
tion, at least using this simple modelling of photoionisation 
feedback. 




Figure 1. The impact of the photo-ionizing background on 
the predicted HI mass function. The panels show models for 
varying Vcut (top panel) and 2cut values (bottom panel). The 
symbols correspond to observations from HIPASS (open trian- 
gles; cf. [Zwaan et al.||2005^ and ALFALFA (filled squares; cf. 
[Martin et al.||2010[ |. We vary the photo-ionization strength be¬ 
tween 25kms“^ ^ Vcut Si 45kms“^ (at a fixed 2cut=10) and 
0 2cut ^ 10 (at a fixed Uut = 30kms“^), as labelled in each 
panel. 


3.2 Other possibilities: cosmology, supernova 
feedback recipe, and star formation law 


In this subsection we explore the effects of a range of addi¬ 
tional aspects of galaxy formation on the HI mass function. 

First, we inspect the impact of cosmology (top panel 
of Fig. [^. The default model has the cosmology used in 
the Millennium N-body simulation of Springel et al. (20051 
as noted in Sec This is compared with a model having 
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the WMAP7 cosmology for which the parameters in GAL- 
FORM have been adjusted to reproduce a suite of other 
observations (Gonzalez-Perez et al. 20141. The abundance 
of galaxies across the full mass range of HI masses in the 
WMAP7 cosmology is larger than in the Millennium cos¬ 
mology. However, the non-monotonic feature described in 
the last subsection remains and therefore is not dependent 
on cosmology. 


Second, we have explored the SNe saturation model sug¬ 
gested by Font et al. (2011 1 . The default GALFORM model 
parameterizes the SNe feedback mass loading efficiency as 
P = (Kirc/Vhot)~“‘'°‘, where Kirc is the circular velocity of 
the galaxy at the half-mass radius. The parameters Idiot and 
Ohot are adjustable and control the strength of SNe feedback. 
The default model has Vhot = 485kms“^ and Qhot=3.2 (cf. 


Bower et al. 20061. In order to ht the number of satellite 


galaxies and their metallicity in the Milky Way, [Font et al.| 
| |2011[ ) proposed a scenario in which the SNe feedback effi¬ 
ciency is not a simple power law for all galaxy masses (as it 
is in the default model). They modelled an efficiency of SNe 
feedback that saturates in low-mass galaxies with circular 
velocities less than 65kms“^, but which follows the default 
model for larger galaxies (see [Lagos et al.|2013 for a phys¬ 
ical interpretation of this saturation). Here we define the 
circular velocity of saturated low-mass galaxies as Idat, and 
explore the impact on the HI mass function of varying 14at 
(from 0 —70kms“^, middle panel of Fig. [^. The abundance 
of galaxies with Mhi < 10 ^ h~^ Mq increases as the value of 
Kat increases, because less gas is reheated and expelled from 
the disk to the halo when the galaxy has circular velocity 
Kirc Si Kat. However, as we found for the case of changing 
cosmology, the non-monotonic feature is still present in the 
predicted HI mass function even after varying Kat. 


In addition, we explore the effect of changing the pa¬ 
rameters of the star formation law within the observed 
range as noted in Sec. The bottom panel of Fig. 
shows the HI mass function resulting from varying val¬ 
ues of izsF = [0.27, 0.5, 0.755] (solid-, dotted-, dashed-lines) 
Gyr“^, Po/fcB=[17000, 30000, 40738] (black-, blue-, red-line) 
cm~^K and /?preas=[0.8, 0.9, 0.99] (the thinest to the thick¬ 
est). We show the HI mass functions from the sets of 27 
combinations of three parameters of the star formation law. 
The varied vsf and Po values lead to an overall shift in the 
HI mass function compared to the default model. Varying 
the value of /3presa also produces an overall shift in the mass 
function. In addition, varying Ppreaa results in a large change 
in abundance at the low mass end of the HI mass function. 
This change is most likely due to a lower T,h 2 /'Shi in the 
regime Pext < Po (HI dominated regions) which is obtained 
for large values of /3press ~0.99. 


Figure 2. Impact of various galaxy formation physics on the HI 
mass function; Top panel: The effect of the changing the cos¬ 
mology from WMAPl to WMAP7. Middle panel: The effect of 
varying the circular velocity below which supernova feedback sat¬ 
urates, as described in jFont et al.|p011| |. Bottom panel: The ef¬ 
fect of varying the parameters in the star formation law from 
[Lagos et al.| ( |2012[ within the observational constraints (see test 
for details). The symbols correspond to data from HIPASS (open 
triangles; cf. jZwaan et al.|2005| l and ALFALFA (filled squares; cf. 
[Martin et al.|2010||. ' ' 


However, as in the cases of varying cosmology and 
SNe feedback, changing the SF law still results in a non¬ 
monotonic HI mass function for Mhi SS 10® h~^ Mq in all 
27 cases, in disagreement with observations. This suggests 
that new or modihed physics is needed in the default galaxy 
formation model in order to improve the agreement with the 
observations at Mhi ^ lO®/i“^M 0 . 
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4 AN IMPROVED MODEL FOR 
PHOTO-IONISATION FEEDBACK 

In the previous section we explored a range of possibilities 
for phenomena which may govern the low and intermediate 
mass range of the HI mass function. In all cases the predic¬ 
tions include a non-monotonic feature which is not present 
in the observations. The results of Fig. suggest that the 
non-monotonic feature in the low mass end of the HI mass 
function is an artifact of adopting a single value of Kut in the 
parametrization of photo-ionisation feedback. We therefore 
model photoionisation feedback using a redshift dependent 
Hcut {z) which is motivated by the calculations of |Sobacchi &l\ 
Mesinger (20131 using ID collapse simulations. Note that in 
our default model, we use /Spress = 0.99 in the star-formation 
law. 


Sobacchi & Mesinger (20131 assumed the following func¬ 
tional form for the critical mass, Merit, defined as the total 
halo mass where the baryon fraction is half of the global 
value: 


Merit — Mq J2I 


1 + Z 
10 


1 - 


1 + z 
1 -h 2 in 


( 3 ) 


where Mo,a,b,c and d are parameters that Sobacchi & 
Mesinger fitted to the results of their numerical simula¬ 
tions. The square bracketed term of Eqn. is related to 
the redshift evolution of the UVB exposure and redshift. 
The parameter zin corresponds to the redshift at which 
the halo was exposed to a UV background. |Sobacchi fc| 
Mesinger (20131 explored a large parameter space of halo 
mass (Mo), UV background intensity (J 21 , expressed in units 
of 10“^^erg s“^ Hz“^cm“^ sr“^), redshift ( 2 ), and redshift of 
UV background exposure (zin) of galaxies. These are uncer¬ 
tain because the details of how reionization proceeded are 
yet to be observed. 


Sobacchi & Mesinger (20131 quoted the best fitting val¬ 
ues of the parameters in Eqn. to their simulation results 
as (Mo, a, b, c, d) = {2.8 x 10® M©, 0.17, -2.1, 2.0, 2.5). For our 
purposes, it is helpful to recast the critical halo mass as a 
circular velocity, Vcut,(Cole et al.|[l994( 


Merit = 3.36 X 10® 


Uer 


kms" 


(l + z) 2 Mo, (4) 


galaxy luminosity function in the optical. We find that the 
model gives a better match to the observed HI mass function 
if we instead use a critical mass at 2 = 0 which is reduced 
by a factor of 10 with respect to the value advocated by 
Sobacchi & Mesinger, corresponding to Kuto ~ 50kms“^. 
This difference could arise due to the ID simulations used by 
Sobacchi & Mesinger missing 3D substructures that favour a 


lower value for the critical mass at z=0 (see Sec 3.1 in Sobac 


chi & Mesinger 20131. In what follows, we treat Muto and 


a„ as free parameters in the new photo-ionisation feedback 
model. 

In the top left panel of Fig. we show the predicted HI 
mass function at 2 = 0 for models using Veut ( 2 ) (Eqn. 
with Veuto = 50kms“^, while varying from -8.2 to -0.2. 
In the top right panel of Fig. we show the predicted HI 
mass function for models using Vmt{z) and q ;„=- 0 . 8 , while 
varying Muto from 35 — 60kms“^ at z=0. The top panels 
of Fig. I^show that this model is able to capture the shape 
of the HIMF down to Mhi ~ 10® h~^ M© much more ac¬ 
curately than the default model does at z=0. This suggests 
that an evolving photoionisation model is preferred by the 
observations in the local Universe. 

We also show the predicted HI mass function at z=0.1 
(middle panels) and z=0.25 (bottom panels) in Fig. In¬ 
terestingly, the predicted HI mass functions at these red- 
shifts show larger variation on changing the parameters of 
Eqn. ([^ than the HI mass function at z=0. Ongoing and 
future Hl-selected galaxy surveys using the SKA (Square 
Kilometre Array; e.g. [Baugh et al.||2004l [Power et al.||201()l 
Kim et al.j2011 1 and its pathfinders, such as ASKAP (Aus¬ 


tralian Square Kilometre Array Pathhnder; cf. [Johnston] 


et al. 20081, and MeerKAT (Meer Karoo Array Telescope; 


cf. Jonas 2007 \ are expected to extend our view of the HI 


Universe to higher redshifts. These observations will probe 
the physical processes that drive galaxy formation, and we 
expect these surveys to be particularly enlightening regard¬ 
ing how photoionisation feedback should be modelled. 

Encouraged by the success of this model in the lo¬ 
cal Universe, we next constrain the parameters of the new 
photo-ionisation model, Muto and For combinations of 
values 6 — (MutOi Q^u)) we obtain the of fh® model with 
respect to the observed HI mass function, and calculate the 
likelihood of each combination of parameters, 


where we have assumed a Hubble constant of 
70kms“^Mpc“^. Combining Eqns. l|^ & (|^, we ob¬ 
tain a redshift dependent Vcut{z), which, after adopting the 
suggested best fitting values for the parameters b, c and d 
for Eqn. is given by 


Kut( 2 :)[kms = Kuto(l + «)“' 


1 + z 
1 -|- 2in 


21 2 . 5/3 


( 5 ) 


Here the Kuto is the circular velocity of dark matter halos 
at z=0 below which photoionisation feedback suppresses gas 
cooling and = 6/3 -1-1/2 in the notation of Sobacchi & 
Mesinger. For the best htting value of 6 = —2.1, = —0.2. 

For simplicity, we assume the redshift of UV background ex¬ 
posure to be ziN=10. The value of Veuto is ~ 100 kms“^ if we 
set J 21 = 1. If we use this value of Kuto, the model under¬ 
predicts the low-mass end of the HI mass function by a sub¬ 
stantial amount and also underpredicts the faint-end of the 


C oc exp > 

, 2 [^i(®) “ Vobs.i]® 

where X = 2 ^- 2 - 

i=\ obs,i 

where Xi{6) and CToUs,* represent the model prediction for 
the i*® observed data point, Xobs.i, and its measurement er¬ 
ror, respectively. From this two-dimensional likelihood dis¬ 
tribution, we estimate the joint 1-a (68.3%), 2-cr (95.4%), 
and 3 -ct (99.7%) conhdence levels around the best combi¬ 
nation of parameters (those that minimize the x^)- This is 
shown in Fig. for and Kuto- 

We find that the combination of parameters, (14uto, au), 
that provides the best fit to the observed HI mass func¬ 
tion at z=0 is (Ucuto,®!)) = (55kms“^, —2.02) for zin=6 
and (50kms“^, —0.82) for zin=10 (hereafter, our best fit 
model). Two values for 2 in are considered to reflect the 
current uncertainty in the epoch of reionization. We find 
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Figure 3. The impact of the redshift dependent Vcut(^) modelling on the predicted HI mass function (top panels) at z=0. The symbols 
correspond to observations in the local Universe from HIPASS (open triangles; cf. |Zwaan et al.|[^005[ | and ALFALFA (filled squares; 
cf. [Martin et aL]|2010[ |. Note that these data are reproduced without error bars in the z=0.1 and z=0.25 middle and bottom panels 
as a reference to illustrate the evolution of the mass function. We vary a.v from -8.2 to -0.2 (which controls the redshift dependence), 
and Ucuto from 35 — 60kms“^, as labelled. The middle and bottom panels show the predicted HI mass functions at z=0.1 and z=0.25, 
respectively, for the same combination of parameters shown in the top panels. 
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Figure 4. The upper panels (a) and (b) show the likelihood contours and the best fitting parameters Qv — Fcuto (indicated by the 
cross) for zij,r=6 and zin=10, respectively. The red, orange, and yellow colours indicate the l-cr, 2-a, and 3 -(t levels. The maximum 
l ikelihood positions are shown as a cross symbols in the figures, (Fcuto, chj) = (55,-2.02) for zjjv=6 and (V(;uto, au) = (50,-0.82) for z/jv=10 
. Sobacchi & Mesinger (2013i suggested a^=-0.2. and Fcuto ~ 55(40)kms“^ (which corresponds to the circular velocity advocated by 
Sobacchi &; Mesigner after reducing their critical halo mass by a factor of ten; see details in fj^ for UV background intensities o f 1 (0.01) 
[10“^^erg/s/Hz/cm^/sr]. The panel (c) shows the ratio of (ns defined in Eqn.[^ between the default model and the model of |Sobacchi| 
I&: Mesinger| | [^13[ l (SM2013) adopting a UVB intensity 0.1, parameters using in the three different confidence levels, and the best set of 
parameters in our modelling for zin=10. The panel (d) shows the evolution of V'cut(2) for |Sobacchi &: Mesinger] | |2013[ l (solid lines, two 
different UV background intensity values [0.01,1]), two examples from the 2-f7 likelihood region (dotted red lines), two examples from 
the 3-0- likelihood region (dashed blue lines), and the maximum likelihood position (long dashed red line) for ziim=10. 


that the best fitting combination of parameters for zin=6 
are in the l-cr confidence level of the two-dimensional like¬ 
lihood distribution for the zin = 10 case, indicating that this 
assumption does not introduce a large uncertainty. These 
values are shown in the panels (a) and (b) of Fig. [^by the 
cross. The likelihoods derived from the modified model cor¬ 
respond to improvements in the predicted low mass end of 
the HI mass function, as demonstrated in the panel (c) of 


Fig. W which shows the ratio of the value (defined in 
Eqn.j^ of the default model (that in |Lagos et al.|20T^ to the 
value of using the Sobacchi & Mesinger model for zin = 10. 
Also shown are the three different representative sets of pa¬ 
rameters drawn from each confidence level in our modelling, 
as well as the best fitting parameter set for zin = 10. The 
panel (d) of Fig. shows the evolution of Vcntiz) for dif¬ 
ferent sets of parameters. The solid black lines show the 
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evolution of Vcut{z) from Sobacchi & Mesinger ([20]^ for 
two different values of the UV background intensit}rj(l and 
0.01 X 10“^^ergs“^ Hz“^ cm“^ sr“^.) and values for Kuto ~ 
of 55 and 40 km s“^. The blue dashed lines in the right panel 
of Fig.|^show the evolution for different combinations of pa¬ 
rameters that lie in the 3 -ct likelihood region, while the red 
dotted-lines show parameter combinations that lie in the 2 -(j 
confidence level. The red long dashed line shows the evolu¬ 
tion for the best fitting set of parameters for zin = 10. 

We note that in this exercise, and 14uto are var¬ 
ied independently of the other parameters in GALFORM. 
We have verified that the modified model for photoionisa¬ 
tion feedback using these parameters does not affect other 
model predictions that were in good agreement with the ob¬ 
servations. This is because these observations tend to probe 
more massive galaxies than the low-mass end of the HI mass 
function. For example, Fig.j^shows the predicted bj and K- 
band luminosity functions in the default model and in the 
new model, demonstrating that we find broad consistency 
with observations for galaxies brighter than -15 mag. 


the abundance of central galaxies drops below the abun¬ 
dance of satellite galaxies in the default model. In contrast, 
our best fit model includes a similar abundance of central 
and satellite galaxies. The reason is that central galaxies 
with circular velocities from 30kms“^ to 50kms“^ in the 
default model are never affected by photoionisation feed¬ 
back, while in the new model these are suppressed between 
0 ^ 2 ^ 51 . These central galaxies dominate the HI mass func¬ 
tion at 10® Mq ^ Mhi ^ lO®h“^M 0 . In addition, cen¬ 
tral galaxies with circular velocities Si 30kms“^ in the de¬ 
fault model cannot contain cold gas at zSilO, and cease 
forming stars shortly after. In the new model these galaxies 
are able to continue accreting new gas and form stars down 
to z~l. This produces the increased abundance of central 
galaxies at Mhi ^ 10^ h ~^ Mq that reside in haloes of mass 
s; 10 ^°/i-^Mq. 

This difference in the abundance of satellite versus cen¬ 
tral HI galaxies explains why the non-monotonic feature seen 
in the default model is not found in our new model, provid¬ 
ing evidence for an evolution in the strength of photoionisa¬ 
tion feedback. 


4.1 The distribution of HI in central and satellite 
galaxies 

In Fig.|^we investigate the contributions of satellite and cen¬ 
tral galaxies to the HI mass function in the default model 
and in our best fit model. The non-monotonic feature seen 
in the HI mass function from the default model is related to 
both the suppression of cooling process in central galaxies, 
and the abundance of satellite galaxies at low HI masses. 
Photoionisation feedback in the default model suppresses 
the cooling of gas in central galaxies which have circu¬ 
lar velocities less than Kut- As a result, the contribution 
of central galaxies to the number density of galaxies with 
Mhi ~ 10^h~^M q is less than 10 % in the default model. 
On the other hand, satellite galaxies in massive dark matter 
haloes (Mhaio > IO^^/i^^Mq) contribute more than 90% of 
the abundance at galaxy masses of Mhi ~ lO®h“^M 0 . Fur¬ 
thermore, in the default model the abundance of satellite 
galaxies monotonically increases towards low HI mass. The 
reason why the HIMF of satellite galaxies does not show the 
dip at low HI masses observed in the HIMF of central galax¬ 
ies is that the satellites were mainly formed before reioniza¬ 
tion. Thus the sum of the abundance of central galaxies and 
satellite galaxies introduces a non-monotonic feature in the 
default model around the Mhi that is directly related to the 
value of the parameter I 4 ut. 

The contribution of satellite galaxies to the HI mass 
function in our best fit model is larger than in the de¬ 
fault model across the whole HI mass range shown. Cen¬ 
tral galaxies dominate the abundance of galaxies at Mhi 
10® h~^ M 0 in both models. However, at Mhi 5; 10^ h~^ M 0 , 


^ [Sobacchi fc Mesinger] ( [2013^ pointed out that their results are 
not very sensitive to the intensity of the UVB: increasing the 
intensity of the UVB by two orders of magnitude increases the 
critical mass only by a factor of 2-3. Measurements of the Ly¬ 
man alpha forest imply that the UVB intensity lies between 
(0.1 — 1) X 10“^^ ergs“^ Hz“^ cm“^sr“^ at z=2 to 6 (e.g. [Bolton] 
[&: Haehnelt[[2007[ [Haardt fc Madau[[T012[ |. But its evolution for 
2 > 6 is very uncertain. 


5 OTHER PROBES OF PHOTOIONISATION 
FEEDBACK 

In this section, we examine the effects of the modified model 
for photoionisation feedback on the predicted spatial distri¬ 
bution of galaxies selected by their HI mass, and on the 
relation between the HI and the stellar mass contents of 
galaxies at z= 0 . 


5.1 Clustering of HI galaxies 


We first use the two-point correlation function to investi¬ 
gate how the modified photoionisation feedback influences 
the spatial distribution of galaxies selected by their HI mass. 
The distribution of Hl-selected galaxies in the default model 


agrees well with the observations of Meyer et al. (20071 and 
Martin et al. (20101 who used the HIPASS and the AL¬ 


FALFA surveys, respectively, to measure the cluste ring of 
galaxies with Mhi > lO®h“^M 0 (Kim et al. 20131. How¬ 


ever, it is challenging to measure the spatial distribution of 
Hl-poor galaxies due to the small number of such galaxies 
observed in current surveys (cf. Papastergis et aL][2013 K 
The top panel of Fig. shows the predicted two-point 
correlation functions measured from the default model and 
for our best fitting model, while the bottom panel shows the 
ratio between the correlation function in the best fit model 
to that in the default model. We consider four HI mass selec¬ 
tions in order to see how the modified model for photoioni¬ 
sation feedback is imprinted on the clustering of Hl-selected 
galaxies. The predicted two-point correlation functions show 
that the clustering amplitude in the default model is much 
lower than in our best fitting model for the HI mass sam¬ 
ples with Mhi > lO®/i“^M0 and Mhi > lO^h“^M0. In 
addition, the slope of the two-point correlation function at 
small separations is steeper in our best fitting model than 
it is in the default model. The ratios of correlation func¬ 
tions for the HI mass samples with Mhi > lO®/i“^M0 and 
Mhi > lO^/i“^M0 shows a large difference at small separa¬ 
tions. This is because the two models have different numbers 
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Figure 5. The bj band and K band galaxy luminosity function at z=0 for the default model (black solid line) and our best fit model 
(blue dashed line). Observations from [Norberg et al.| ( [?002| | for bj band luminosity function (left panel) and [Cole et al.| | |2001[ | and 
[Kochanek et al!] < |2001[ | for K band luminosity functions (right panel) are shown by symbols with errorbars. 
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Figure 6. The HI mass functions at z=0 (black solid lines) in the default model (left panel) and the best fit in our new model (right 
panel). The contributions to the HI mass function, from central galaxies is shown by blue dotted lines and from satellite galaxies as red 
dashed lines. 


of satellite galaxies with HI masses in the range 10® 

and which dominate the small scale clustering 

signal. 

The increase in the observed number of galaxies with 
Mhi < expected from ongoing and future HI 

galaxy surveys (eg. ASKAP, MeerKAT, SKA) promises to 


allow the nature of photoionisation feedback to be inferred 
from the spatial distribution of galaxies with low HI mass. 

5.2 Stellar mass — HI mass relation 

Fig .[^ shows predictions for the relation between the HI and 
stellar mass contents of galax;ies for both the default model 
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Figure 7. The impact of the modified photoionization feedback 
model on the two-point correlation function. The lines show the 
predictions for four galaxy samples defined by different HI mass 
thresholds (see key). The default model (our best fitting model) 
is shown by the thin red (thick black) lines. The lower subpanel 
shows the ratio of the prediction from the best fitting model to 
that from the default model for different HI mass thresholds. 


and our best fit model. These are compared with the ob¬ 
servational results for the median HI mass as a function of 
stellar mass presented by Maddox et al. (20151 (filled cir¬ 
cles). Note that we only consider model galaxies which have 
a HI mass > a stellar mass > and a 

bulge-to-total stellar mass ratio < 0.5 in this plot, in order 
to approximately imitate the selection criteria of the AL¬ 
FALFA (N. Maddox, priv. communication). As argued by 


Maddox et al. (20151, this relation provides a fundamental 


benchmark for galaxy formation models by combining HI 
mass data with stellar mass data. The relation between the 
HI mass and stellar mass in both models remains in reason¬ 
able agreement with the l-cr dispersion of the observations. 
We note that high stellar mass galaxies in the models typi¬ 
cally lie slightly below the observations. This could be due 
in part to remaining differences between the observational 
sample selection and the crude cuts applied to the model 
galaxies to mimic that described above. Indeed, HI observa¬ 
tions of stellar-mass selected samples show somewhat differ¬ 
ent behaviour and are in better agreement with the model 
predictions (Catinella et al. 2010 Serra et al. 20121. Fur¬ 
ther analysis of the HI gas fraction distribution function of 
galaxies at z~0 from GALFORM and comparison with ob¬ 
servations are described in Lagos et al. (20141. 

Finally, we note that the model predictions shown in 
Fig .[^are the true stellar masses rather than those estimated 
from photometry or spectra. Mitchell et al. (20131 showed 
that the process of estimating the stellar mass of a galaxy 
can lead to differences from the stellar mass predicted di¬ 
rectly by the model. These include systematic shifts, which 
arise from a mismatch between the choice of stellar initial 
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Figure 8. The HI mass as a function of the stellar mass for the 
default model and our best fit model with only galaxies which 
have HI mass > stellar mass > and the 

bulge-to-total stellar mass ratio < 0.5. The median HI masses in 
bins of stellar mass are shown for the default model as red empty 
triangles and our best fit model as black empty squares. The 
shaded regions represent the 16% and 84% range for the default 
model (red dashed) and our best fit model ( black solid). The filled 
circles represent the observations in [Maddox et ar] | |2015| |. Blue 
points are derived only from galaxies with SDSS spectra, and 
green points are derived only from galaxies without spectrum. 
The error bars correspond to the l-cr dispersion of observation. 
In comparison, the median HI masses in bins of stellar mass are 
shown for the default model as red filled triangles and our best 
fit model as black filled squares using lower HI mass threshold 
(Mhi> 1O®/i-2M0). 


mass function used to estimate the stellar mass and that 
adopted in the semi-analytic model and from differences in 
the modelling of dust extinction, and scatter, due to varia¬ 
tion in the assumed metallicity and star formation history. If 
we replaced the true model stellar masses in Fig. with esti¬ 
mated stellar masses, this would soften the downturn in the 
relation at high stellar masses. We also include the predic¬ 
tions from the two models using a lower HI mass threshold 
(Mhi>10®/i“^Mq instead of Mhi>1O®/i~^M0; see red filled 
triangles for the default model and the black filled squares 
for our best fitting new model connected by long dashed lines 
for the median HI masses in bins of stellar mass). The pre¬ 
dictions deviate from one another at stellar masses between 
lO^/i“^M 0 and lO®/i“^M 0 . This stellar mass range corre¬ 
sponds to HI masses between lO^/i^^M© and lO^/i^^M©, 
where the two models also predict different HI abundances 
because of the different modelling of photoionisation feed¬ 
back. 
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Figure 9. The HI density in units of the critical density as a 
function of redshift for our Vcut ( 2 ) model (this work, black solid 
line) and |Lagos et al.] l |2012[ | model (red dashed line). The points 
show various observations. The parenthesis after each observa¬ 
tion indicated the method used are HI: direct 21cm line detection 
I Zwaan et al.|[2005[ [Martin et al.||2010[ [Freudling et al.||2011[|, 
St.: stacking of 21cni line [Lah et al.|2007[ [Delhaize et al'|2013[ 

Rhee et al.[[2013[, D LA: absorption line stu dies ([Peroux et al.[ 
2003[[Rao et al.|2006[ [Noterdaeme et al.|2012[ [Zafar et al.|2013[l, 

and i.m.: intensity mapping of the 21cm line ( [Chang et al. ]20 To| 

[Masui et al.|2013[ |. 


5.3 Evolution of HI density 


Finally we investigate the impact of our new model on the 
HI density as a function of redshift. Fig. shows Hhi, de- 
Hned as the ratio between the mean HI density and the crit¬ 
ical density of the Universe, for the new evolving 14ut(2) 


model (black solid line) and the Lagos et al. (20121 model 
(red dashed line, the default model). The HI density pre¬ 
dicted using the new model with I4ut(z:) shows slightly bet¬ 
ter agreement with the observed HI density at z<0.1 than 


the Lagos et al. (20121 model, recovering the decreases Hhi 


suggested by the observations of z^O.2. However, the differ¬ 
ence between the two models are within the uncertainties of 
the observations, and new constraints at intermediate red- 
shifts are needed. 


6 SUMMARY 


One of the most important ingredients in a model of galaxy 
formation is the description of neutral hydrogen which pro¬ 


vides the fuel for star formation. The HIPASS (Zwaan et al. 
|2005| ) and ALFALFA ( [Martin et al.|20lo| ) surveys have been 
used to estimate the HI mass function down to masses of a 
few times 10® h~^ Mq in the local Universe. Attempting to 
reproduce these observed HI mass functions using a galaxy 
formation model allows us to probe the physics of low mass 
galaxy formation. 

In the standard implementation of GALFORM the ef¬ 
fect of photoionisation feedback is modelled by restricting 
galaxy formation in halos which correspond to effective cir¬ 


cular velocities smaller that a Hxed circular velocity of 14ut 
at redshifts following the end of reionization Zcut (the default 
model). Halos with circular velocities below 14ut are not al¬ 
lowed to cool gas at redshifts below Zcut- The predicted HI 
mass function from this default model shows good agreement 
with observations of the HI mass function for masses greater 
than Mhi ~ 10® h~^ Mq. However, the default model fails to 
explain the abundance of galaxies with HI masses between 
1 O®/i- 2M0 and 1 O®/i- 2M0. 

There are several possibilities to ease the discrepancy 
at the low mass end of the HI mass function between ob¬ 
servations and the default model, including chaning the su¬ 
pernova feedback, cosmological parameters and the form of 
the star formation law in the model. In particular [Kim et al.| 
(20131 showed that photoionisation feedback is an important 
physical processes for the low mass end of HI mass function. 
We have investigated these possibilities using a version of 
the GALFORM semi-analytical galaxy formation model (cf. 
Gole et ar]|2000 1 by Lagos et al. (20121 to simulate the low 


mass end of the HI mass function. In our investigation we 
showed that varying the SNe feedback, the efficiency of star 
formation and the cosmological model all failed to explain 
the shape and abundance of the low mass end of HI mass 
function. Varying the redshift and critical velocity of galaxy 
suppression of formation by photoionisation feedback from 
reionization is also unable to explain the shape of the low 
mass end of HI mass function. 

To understand what photoionisation physics affects the 
shape and abundance of the HI mass function in the local 
Universe, we introduced a redshift dependent photoionisa¬ 
tion feedback into GALFORM, motivated by the simula¬ 
tions of Sobacchi & Mesinger (20131. This illustrates the 
modular nature of semi-analytical modellling, whereby the 
treatment of a given physical process can be overhauled and 
improved when new information becomes available. We find 
that redshift dependent modelling of feedback from pho¬ 
toionisation on low mass galaxy formation is needed to ex¬ 
plain the shape and abundance of the observed HI mass 
function in the local Universe. We also find that the sensi¬ 
tivity of the HI mass function to the redshift evolution of 
photoionisation feedback strength is larger at high redshift. 
Our modelling suggests that future measurements of HI clus¬ 
tering in low mass galaxies and the relation between the HI 
mass and the stellar mass of galaxies will provide additional 
constraints on the form of ionising feedback. 

Ongoing and future HI selected galaxy surveys on the 
SKA and its pathhnders, such as ASKAR, MeerKAT are 
expected to extend our view of the HI Universe to higher 
redshifts not only for the HI mass function but also for the 
distribution of HI galaxies. The evolution of the HI mass 
function as a function of redshift will tell us about the phys¬ 
ical processes that drive low mass galaxy formation, includ¬ 
ing photoionisation feedback. 
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